Quintessence and the Relic Density of Neutralinos 



a,b * 



Pierre Salati 

a) Laboratoire de Physique Theorique LAPTH, B.P. 110, F-74941 Annecy-le-Vieux Cedex, France. 

b) Universite de Savoie, B.P. 1104, F-73011 Chambery Cedex, France. 



July 17, 2003 - LAPTH-951/02 



The archetypal model for the recently discovered dark energy component of the universe is based 
on the existence of a scalar field whose dynamical evolution comes down today to a non-vanishing 
cosmological constant. In the past - before big-bang nucleosynthesis for that matter - that scalar 
field could have gone through a period of kination during which the universe has expanded at a 
much higher pace than what is currently postulated in the standard radiation dominated cosmology. 
I examine here the consequences of such a period of kination on the relic abundance of neutralinos 
and find that the latter could be much higher - by three orders of magnitude - than what is estimated 
in the canonical derivation. I shortly discuss the implications of this scenario for the dark matter 
candidates and their astrophysical signatures. 



I. INTRODUCTION. 

The recent WMAP observations of the Cosmic Microwave Background (CMB) anisotropics [1], combined cither 
with the determination of the relation between the distance of luminosity and the redshift of supernovae SNela [2] , 
or with the large scale structure (LSS) information from galaxy and cluster surveys [3], give independent evidence for 
a cosmological average matter density of Om = 0.27 ± 0.04 [1]. This value may be compared to a baryon density of 
Ob = 0.044 ± 0.004 as indicated by nucleosynthesis [4] and the relative heights of the first acoustic peaks in the CMB 
data. A significant fraction of the matter in the universe is dark and non-baryonic. The cosmological observations 
also point towards a flat universe with Otot = 1-02 ± 0.02 and strongly favour the existence of a cosmological constant 
which contributes a fraction Oa = 0.73 ± 0.04 to the closure density. The pressure-to-density ratio w of that fluid is 
negative with a value of w = — 1 in the case of an exact cosmological constant. That fl\ component is called dark 
energy as opposed to the Om dark matter contribution. 

The nature of the astronomical dark matter is still unresolved insofar. The favorite candidate for the non-baryonic 
component is a weakly-interacting massive particle (WIMP). The so-called neutralino naturally arises in the frame- 
work of supersymmetric theories. Large efforts have been devoted in the past decade to pin down these evading 
species. New experimental techniques have been devised to look for the direct and indirect astrophysical signatures 
of the presence of neutralinos in our Milky- Way [5] as well as in extra-galactic systems [6]. The uncertainty on the 
theoretical estimates of the various signals has been considerably reduced. As an illustration, the energy spectrum of 
secondary spallation antiprotons - the natural background to a putative neutralino-induced antiproton extra radia- 
tion - is now well under control [7] . Another example of the level of sophistication which the theoretical investigations 
have reached is provided by the calculations of the neutralino relic abundance O^^. The observation that this relic 
density - depending on the numerous parameters of the model - falls in the ballpark of the measured value for Om 
has been a crucial argument in favor of supersymmetric particles as a viable option to non-baryonic dark matter. 
A large number of processes - typically ~ 2000 - are now taken into account and the corresponding diagrams are 
automatically generated and calculated with the help of numerical codes such as micrOMEGAs [8] . Co-annihilations 
are taken into account and the thermal averaging (av) of the product of the velocity by the cross section is performed. 

Surprisingly enough, calculations of fl^ are based on the assumption that the universe is dominated by radiation 
when neutralinos decouple from the primordial plasma and reach their relic density. This hypothesis is presumably 
correct as soon as primordial nucleosynthesis (BBN) sets in at a time of ~ 1 second. We have however little information 
on the earlier pre-BBN period, a crucial stage during which neutralinos freeze out. If the expansion of the universe 
is modified with respect to a pure radiation-dominated behaviour, the quenching of these species could be drastically 
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modified. An increase in the expansion rate H accelerates the decouphng of neutralinos and translates into larger 
values for the relic density Jl^. 

Exploring the effects of a modified expansion history of the universe onto the relic abundance of neutralinos is 
no longer a mere academic exercise. Such an analysis has become mandatory inasmuch as a new and unexpected 
component - the dark energy J^a ~ has been discovered. The potential interplay between that component and its 
matter counterpart f2m is worth being explored and may have unexpected consequences. The archetypal model for 
the cosmological dark energy is the so-called quintessence [9,10] and relies on the existence of a neutral scalar field 
<t> with Lagrangian density 

£ = ^9^''d^'Pd,.<P - . (1) 

Should the field <!> be homogeneous and the metric be flat, the energy density may be expressed as 

p^^n = ^ + , (2) 

whereas the pressure obtains from T!y = —g ij P so that 

P$ = ^ - V{^) . (3) 

If the kinetic term <I'^/2 is negligible with respect to the contribution of the potential V"(<f>), a pure cosmological 
constant with w$ = = —1 is recovered since p<i, = — P$ = V ($). As indicated by cosmological observations, 

this is the case today. But the field $ has been continuously rolling down. Should the kinetic term <i>^/2 have 
dominated over the potential V (<&) in the early universe, a period of kination - i.e., domination by the kinetic energy 
of the field $ - would have ensued with drastic effects on the expansion rate of the universe [11]. 

In section II, we briefly recall why a pure cosmological constant should be disregarded and replaced by a dynamical 
dark energy component in the form of a scalar fleld, the so-called quintessence whose salient features are presented. 
The existence of tracking solutions provides a natural solution to the problem of initial conditions. We also pay 
some attention to the difficulty of generating a kination-dominated expansion in the early universe together with 
a cosmological constant today [10]. We show that this difficulty may be circumvented depending on the potential 
V ($) that drives the evolution of the scalar field and we propose examples where quintessence boosts the expansion 
rate in the past while it still accounts for today. Following a suggestion by [12], we investigate in section III the 
effects of kination on the thermal decoupling of neutralinos and derive an approximate relation between their relic 
abundance fi^ and their annihilation cross section in the presence of kination. Section IV is devoted to a discussion 
of the consequences of this scenario on the astrophysical signatures of neutralino dark matter. 



II. KINATION AND THE EXPANSION RATE OF THE UNIVERSE. 



Two difficulties arise with a pure cosmological constant. The coincidence or fine-tuning problem lies in the fact 
that the vacuum energy comes into play only today and is therefore of order the closure density p^, an exceedingly 
small value with respect to the typical Planck energy scale -/V/pjg^j^j.]^. set by particle physics 

Pa . 4keVcm-3 . ^ (4) 



M^lanck {1.22 X 1019 GcV}" 

The other issue is related to the initial conditions in which a pure cosmological constant has to be prepared. At the 
Planck time, the corresponding vacuum energy density p\ = p^ needs to be exceedingly fine-tuned with respect to 
the radiation density 



Pjad P?ad UfrlanckJ 



10-''" . (5) 



Because a pure cosmological constant does not vary in time, similar values are obtained in the previous relations. 
This has led to some confusion between the fine-tuning and the initial condition problems. Quintessence actually 
solves only the latter difficulty while the fine-tuning of p° with respect to Mp^^^^^ is forced by direct observation. 
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The idea of quintessence is based on the existence of a scalar field $ that rolls down its potential V according to 
the homogeneous Klein-Gordon equation 

dV 

$ + 3i7$ + - . (6) 

Should $ satisfy to that relation, its energy-momentum tensor would be immediately conserved as may be inferred 
from the identity 



D^T^" - (a» • {D^id^^) + ■ (7) 



That energy-momentum tensor obtains from the Lagrangian density (1) and may be written as 

Tf^„ = df^ipdyip - Qf^iyC . (8) 
Its conservation translates into the adiabatic expansion of quintessence as 

where the dark energy density p$ and pressure P$ have been defined in relations (2) and (3). The expansion rate 
H = d/a is increased by the presence of quintessence 

= ^ {PB + P$} , (10) 

where the energy density of the background pB = Piad + Pm is dominated by radiation in the past and by matter since 
equality. Notice that at fixed equation-of-state w — P/ p, the energy density scales as 

poca-3(l + w) . (11) 

As discussed in [10], we are looking for tracking solutions of (6) for which wb > ~ constant. This allows a natural 
and smooth decrease of the quintessence energy density p$ that remains subdominant until recently when it takes over 
from the background ps- This scenario is realized only if IF'/^I decreases with V as time goes on. This translates 
into the condition 

V"V 

r = -^>i. (12) 

where the prime denotes the derivative with respect to the field As long as the background dominates - for 
ri$ — /9$/(p$ + pb) ^ 1 - the parameter T is related to the equation-of-state of the tracking solution through 

- 2 (r - 1) 
"'^ - i + 2(r-i) ■ 

Inverse power law potentials V oc <I>^" with a > naturally drive quintessence since F — 1 = 1/a > 0. On the 
contrary, exponential potentials V cx exp (— $/M) for which F = 1 and w$ = wb should be disregarded. 

As shown in [10], the tracking solutions are attractors towards which the field $ relaxes for a large range of values of 
the initial energy density p]^. Quintessence solves therefore the initial condition problem to a large extent - actually 
to the extent that p^ is comprised between p° and pg . The Klein-Gordon equation (6) may be written as 



V 



1 dine 



= 3,/-VTT^jl + -^^ , (14) 



where k = 8 tt G / 3 and 

$2/2 1 + 
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Along the tracking solution, and x are constant. When quintessence pops above the background - when p$ pB 
- relation (14) translates into 



A-Zpianck V 24 TT V 1 + 

in the case of an inverse power law potential. The field $ is naturally of order the Planck mass Mpianck when 
quintessence takes over from the background. The actual coincidence - forced by direct observation - lies in the fact 
that quintessence becomes dominant today. As its energy density - and therefore the potential V ~ are of order the 
closure density p^, this implies a significant amount of fine-tuning as 

V{^o) = / lO^^^GeV^ , (17) 

where $o ^ -^^Pianck is the present value of the scalar field. We conclude at this stage that quintessence does not solve 
the coincidence or fine-tuning problem [13]. This point cannot be too strongly emphasized. 

As also sketched in [10], the entire scenario of quintessence would be in jeopardy should the initial energy density 
p]j, exceed pg. If so, the initial position from which the scalar field starts to roll lies too high on the potential. The 
field $ falls precipitously. As its kinetic energy becomes dominant, a kination stage ensues during which the potential 
V plays no role. The Hubble friction - see the Klein-Gordon equation (6) - slows down the fall and the field gets 
frozen at the value <i>F such that 

T7^^yj(l + |lnaj , (18) 
Aipianck V 4 7r ( 2 ) 

where f^i = p'g./p'g- If, in the early universe, quintessence drives a period of kination with fii ;:!> 1, the field gets frozen 
at $F ~ 10 — 100 Mpianck, & valuc that exceeds by far what is required by the conventional quintessence scheme 
where $o ^ A^Pianck- The field has overshot the tracking solution and its energy density has become so small - the 
potential V decreases with increasing $ - that it cannot account for p^. A period of early kination seems therefore 
incompatible with a near constant dark energy density that becomes dominant today. 

The problem raised by overshooting lies in the incompatibility between the conditions $f S> A/pianck and $0 ^ 
Afpianck- However the latter has been derived in the framework of inverse power law potentials. We anticipate that 
the overshooting difficulty could be circumvented with other potentials. As an illustration, we consider the toy-model 
where 

y(<i>) = / exp/^l , (19) 



with M ^ O (A/pianck)- In order for dark energy to overcome the background only today, we still need - as before - a 
considerable amount of fine-tuning as /i^ ~ ~ 10~^^ GeV^. That potential has nevertheless the correct behaviour 
with r = V"V/V''^ = 1 + 2 $/M > 1 and an associated equation-of-state 

= - {1 + ' ' (20) 

on the tracking solution. The initial condition problem is solved to the extent that p]^ is now larger than pg and 
consequently drives a period of kination. In that case, the field rapidly freezes at ^ -Mpianck ~ M - irrespective 
of its precise initial value - and the potential reaches an asymptotic value of 

V ($F » M) ~ p4 - pV (21) 

This example is remarkable as early kination is now mandatory in order to get today a cosmological constant with 
w$ ~ — 1. 
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In our counter-example to the overshooting no go theorem, the potential (19) has a sharp decrease for small values 
of the field and then varies smoothly for $ 3> Afpianck- A dynamical realization of that idea - proposed in [14] - is 
based on 

$2} = M"+4 {$i<i>2}-"/2 , (22) 

where two scalar fields come into play. If one of the fields starts not too far from the tracking solution whereas the 
other one is dropped from an elevated position on the potential, we still get kination in the past whereas a cosmological 
constant is recovered today. This behaviour occurs for a wide range of initial conditions on $1 and <I>2. Other examples 
are presented in [14] with the same trend and we conclude at this stage that a period of early dominant kination could 
perfectly occur at the time of neutralino decoupling without precluding a significant contribution today from a near 
constant dark energy density. 

In the standard big bang model, the early universe is only filled with a gas of ultra-relativistic particles that behaves 
like a radiation of photons. If the kinetic energy $^/2 of some additional scalar field component dominates both the 
potential V ($) and the radiation energy density prad, a period of kination sets in. The overall energy density decreases 
consequently like 

$2 _ 

Ptot ~ p$ ~ — cx a , (23) 

with respect to the scale factor a. This amounts to say that the derivative dV/d^ is negligible with respect to the 
damping term 3 H $ in the Klein-Gordon equation (6) as the field tumbles down the potential. The energy density 
of radiation is related to the temperature T of the primordial plasma through 

Prad = 5cff(r)^T4 _ (24) 

The effective number ^cff of degrees of freedom allows to express piad in units of the corresponding photon density p^. 
The evolution of the scale factor a with respect to the temperature T follows from the requirement that the entropy 
of the radiation is conserved so that 

<x ^Kfi{T)T^a^ (25) 
remains constant as T drops. The effective number hcs of entropic degrees of freedom is very close to gcs since 

/.c.c.5c...E?+E^ . (26) 

B F 

The sum is carried over the bosonic gB and fcrmionic gp spin states that correspond to those species that are massless 
at temperature T. We infer that at any given time, the scale factor a and the temperature T are related through 



^cff (To) 1 ^^'^ To 



(27) 



where ao is the scale factor at some temperature of reference To. During the decoupling of ncutralinos, g^s and hes 
do not change drastically so that a varies approximately like T~^ at that time. This will prove to be helpful when we 
derive an approximate formula for the relic density fi^ in the next section. Numerical results are nevertheless based 
on equation (27). 

We parametrize the contribution of quintessence to the overall energy density through 

^* = 4 ■ (28) 

The ratio 77$ of the quintessence-to-photon energy densities is defined at temperature To. The latter has been set 
equal to 1 MeV. Because quintessence should not upset the conventional BBN scenario, we expect in principle r7<j < 1. 
Note however that a non-vanishing value for rj^ may still alter the pre-BBN era because quintessence may be the 
dominant form of energy at that time if kination holds. Its energy density varies actually like 

f feeff(T) )^ fry 
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whereas for radiation, the temperature dependence is smoother with 

Going back in time, rises more steeply than prad- Even for small values of ry$, quintessence dominates over the 
radiation at early times. Up to a numerical constant, we conclude that the expansion rate H is increased by a factor 
of y/rf^ {T/Tq) with respect to the conventional radiation dominated cosmology. The next section is devoted to the 
consequences of such a change in H on the decoupling of neutralinos. 



III. THE FREEZE-OUT OF NEUTRALINOS. 



At early times - as long as the temperature T exceeds their mass rriy^ ~ neutralinos are in chemical equilibrium. They 
steadily annihilate into lighter species while the reverse process is concomitantly active. The annihilation-production 
reaction 

X + X ^ n ,W+W' , Z"" Z\ HH ... (31) 
reaches its thermodynamical equilibrium. This implies a neutralino density 



n 



X 



9xT-^e-y (32) 



_ „ 7.3 -y ( JL 



that only depends on the mass-to-temperature ratio y = vn^jT as well as on the number g-^ of spin states. As soon as 
T drops below as the result of the overall adiabatic expansion, the neutralino population becomes non-relativistic 
and the annihilations take over the thermal productions. Neutralinos are severely depleted until their density n° is 
so low that they fail to annihilate with each other. Reaction (31) stops to be in equilibrium. The resulting freeze- 
out occurs typically for values of the mass-to-temperature of yp ~ 20. Because the probability for a neutralino to 
encounter a partner has become less than unity per typical expansion time the density n^^ remains subsequently 

constant per covolume - per volume that expands with the expanding universe. The relic abundance fi-^ readily 
obtains from the value of at freeze-out. 

Assuming that there are as many species x than antiparticles x ~ which is an obvious statement if neutralinos are 
Majorana fermions as is the case for instance in the supersymmetric extensions of the standard model - the density 
evolves according to the basic equation 

^ = -3Hn^- < CTanW > ?^x^ + < '^^nV > ■ (33) 

The first expression in the right-hand side refers to the dilution resulting from the expansion of the universe. The 
second term accounts for the neutralino annihilations. The last expression - for which detailed balance has been 
assumed - describes the back-creations of XX pairs from lighter species. The neutralino density is given by as 
long as a thermodynamical equilibrium is reached for reaction (31). In the non-relativistic regime at stake here, it is 
given by relation (32). In terms of the codensity = Uy^/T^ , the evolution equation simplifies into 



dfx 



{< (Tanl- > n^} /x = < fTanf > /° ' . (34) 



As discussed in the previous section, the temperature T varies roughly as the inverse of the scale factor a. The 
codensity corresponds therefore to the number of particles inside a typical volume oc that follows the expansion 
of the universe. In order to solve equation (34), two typical time scales may be defined. To commence, as long as 
reaction (31) has reached its equilibrium, the time derivative df^/dt may be neglected and we recover = as 
we should. The characteristic time scale with which the neutralino codensity relaxes towards its kinetic - and 
thermodynamical - equilibrium value is set by the annihilation rate 

= < ^an" > ""x ■ (35) 

Then, the time scale of the variations of the equilibrium itself may be expressed as 



eq 



-|L„g{/rr3}.2i/{!^i} , (36) 
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in the non-relativistic regime. As is clear in Fig. 1, the neutraUno freeze-out proceeds in two stages. 

(i) At high temperature, as long as Trei ^ T^q, has plenty of time to relax towards the equilibrium value 
that evolves at a much slower pace during this stage. As a consequence, the annihilation reaction (31) reaches 
thermodynamical equilibrium. A very good approximation for the neutralino density is provided by as featured in 
Fig. 1 where the various neutralino density curves all merge for y < 10 whatever the value of the kination parameter 
77$. As illustrated by Fig. 2 and 3, relaxation becomes progressively less efficient as long as the temperature decreases. 
The freeze-out takes place at the very moment when Trei crosses Teq. 

(ii) From that moment on - below the freeze-out temperature Tp - the relaxation time Trei well exceeds r^q. Whilst 
drops down and vanishes, still decreases a little bit and reaches an asymptotic value of f^^'^ that leads directly 

to the relic abundance il^ h? . 




10° 10^ 10^ 10^ 

FIG. 1. Neutralino codensity as a function of the mass-to-temperature ratio y = m^/T for three different values of the 
kination parameter. For = 0, we recover the standard radiation dominated cosmology whereas for tjs = 0.01 and 7?$ = 1, 
the expansion rate H is significantly increased. This leads to an earlier decoupling and to a much larger asymptotic value for 
the neutralino codensity. 

To illustrate our discussion, we have considered a generic bino-like species with mass = 250 GeV and tan/3 = 50. 
The thermally averaged annihilation cross section has been approximated by 

< (TanW > — a + hx , (37) 

where the parameter x stands for the temperature-to-mass ratio T/m^. We have chosen the values a = 0.1 pb and 
& = 21 pb for the annihilation cross section. Before freeze-out, the neutralino density n-^ is set equal to n^. As soon 
as decoupling has taken place, we integrate the Boltzmann equation (33) in terms of the codensity F^. = /^/'^ where 
the coefficient k.{T) — hcff{T)/hctf{TQ) accounts for the reheating of the radiation as its massive species annihilate 
and are converted into lighter populations. As is clear from equation (27), the codensity evolves exactly as a~^. 
For convenience, we have plotted f-^^ instead of as a function of the m^/T ratio in Fig. 1 for three different values 
of 77$. In Figs. 2 and 3, the age of the universe tu as well as the typical time scales Trei and Teq are also featured as 
a function of the mass-to-temperature ratio m^/T . The kination parameter 77$ has been respectively set equal to 
and 1. When the scalar field dominates over the radiation - see Fig. 3 - the expansion of the universe is accelerated 
with much smaller values for tu- The evolution of Trei with y is not affected by kination. On the contrary, Teq is 
much smaller than in the conventional radiation dominated cosmology. Freeze-out is therefore reached at a higher 
temperature with yp = 13.5 for 77$ = 1 instead of yp = 22.7 when = 0. The decoupling temperature increases 
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FIG. 2. The age of the universe tu as weU as the typical time scales Trei and r^q are featured as a function of the 
mass-to-temperature ratio y = m^/T. The freeze-out occurs at yp = 22.7 when Trei overcomes T^q. The standard radi- 
ation dominated cosmology is assumed here with rj,s, — Q so that tu evolves like y^ . 
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FIG. 3. Same as in Fig. 2 with a kination parameter of = 1. The scalar field dominates over the radiation and the 
expansion is accelerated with respect to the conventional situation. This implies smaller values for the age tu that evolves now 
like y^. The freeze-out point is therefore reached earlier with yp = 13.5, hence a larger neutralino relic density. 

by a factor of ~ 2. This mild variation nevertheless implies a strong rise in the neutralino relic density as a result of 
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the strong dependence of on y - see relation (32) . Fig. 1 features a growth of the asymptotic codensity /^"^ by a 
factor of ~ 350 when 77$ is varied from to 0.01. When 77$ is set equal to 1, that increase reaches a factor of ~ 2900. 

The previous example illustrates the large increase of the neutralino relic abundance when kination takes over 
radiation in the pre-BBN period. The effect of 77$ on Uy, may also be understood in the framework of the simple 
approximation which we derive now. Integrating the cross section (Tan^ over the thermal distribution of neutralinos 
is beyond the scope of this work. We will simply assume here that relation (37) applies. This is certainly correct as 
long as the annihilation does not proceed through a s-channel resonance for which the presence of a pole implies an 
unusual behaviour of the cross section and we will postpone this problem to a later analysis. The freeze-out occurs 
at temperature Tp when Trei ~ T(,q. This condition translates into 

\h + hxF]TF'' If = ^ypHiTp) . (38) 



At temperature T, the expansion rate H is increased with respect to the conventional situation by a factor of Vl + ax^ 
where x = T/m^. The parameter a describes the relative contribution of quintessence to the overall energy density 
as compared to the radiation 



5cff(r) t To 

The evolution of after decoupling is followed by integrating relation (34) from x ~ xf down to a; = while 
neglecting the right-hand side term. This leads to a decrease of the codensity just after freeze-out from Jf down 
to its asymptotic value with 

f^^ = — . (40) 

The decrease factor fi depends on xf and on the parameter a which we evaluate at the decoupling temperature Tp. 
This assumption is reasonable because most of the evolution of takes place immcdiatly after freeze-out as featured 
in Fig. 1. We readily obtain 

2 , i aA + bxF b] , , 

Ai(xF,a) - 1 = — VI + aaiFM = ) ' (41) 

XF \ ~a + bxF J 

where the factors A and B depend on the combination u = ^/oixf through the relations 

In \u + + 

A{u) = ^ (42) 

u 

and 

5(^) - ^^±4-1 . (43) 

In the conventional scenario for which a = 0, the previous expressions simplify into A{0) ~ 1 and B{0) ~ 1/2. The 
freeze-out codensity may be derived from the condition (38) 



where Alp denotes the Planck mass. The asymptotic codensity f^^^ is defined by relation (40) and we infer for the 
present epoch a neutralino relic density of 

The factor (4/11) (1/k) accounts for the reheating of the photon background that occurs in the period extending 
from the neutralino freeze-out until now. Most of the relevant dark matter candidates have a mass in the range 
between 100 GeV and 10 TeV. Notice also that the freeze-out temperature-to-mass ratio xf ~ 0.1 — 0.2 is not very 
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sensitive to the quintessence parameter ry^. We can safely take gcsiTp) ~ hcs{Tp) ~ 45 for a numerical estimate of 
the previous expression in order to get 



n«ni w _3 yrr^^ r 3 X 10-27 cm3 s-n 

Py ~ 0.69 kcV cm ^ < } . (46) 

p.XF^ I <Cra„i;> j 

That result is to be compared to the closure density 

= 1.879 X 10"^^ g cm"3 /i^ ~ 10.6 keVcm'^ /i^ , (47) 
hence a neutralino relic abundance that may be expressed as 

n^h^ 6.6 X 10- rax^o:^^^! 

where the annihilation cross section is taken at decoupling and may differ from its present value if b well exceeds a. 
The approximation (48) tends to overestimate the relic abundance but is still acceptable. In our generic example 
illustrated in the previous figures, we have numerically derived a value of fi^ h'^ = 0.145 for 77$ = and of il^ = 440 
for ?y$ — 1 whereas expression (48) respectively gives il.^ h? = 0.16 and 465. 



IV. DISCUSSION AND PROSPECTS. 



We have shown that the neutralino relic abundance increases if a period of kination takes place during the freeze- 
out of the species. We derive now an estimate of the corresponding boost factor as a function of 77$ and m^. If a 
dominates over bxp vci the expression of the annihilation cross section, we may even simplify further relation (48) in 
order to get 

^ r2 X 10-27 cm3 s-M , , 

^xh^ - -~ (49) 



in the conventional radiation dominated cosmology. We have taken a value of yp ~ 20 for the mass-to-decoupling 
temperature ratio in that case. On the contrary, if quintessence is the dominant form of energy with a large value for 
the neutralino fossile abundance becomes 



^ ln(2u) \ h ' ^ ' 



where u — ^/axp has already been defined. With a value of yp ^ W - see the previous section - this implies 

u ~ 1.5 X lO"' ^/rj^ mioo , (51) 

so that the logarithm yields a contribution ^ 10. The parameter toiqo denotes the neutralino mass in units of 100 
GcV. A crude estimate of the relic abundance in this regime ensues 

,2 



fix ^ \/rj^ mioo I T 1 . (52) 

We derive a boost factor of ~ 10'^ y/rj^ n^ioo with respect to the conventional cosmology. If now b xp is the leading 
term as regards the annihilation cross section, we find that the relic abundance which is normally given by 

n,,' . . (53, 

is increased to 

o 1.2 ^ r 1.2 X 10-22 cm^s-n , 

il^ ~ y/rj^ mioo < ~ > ■ (54) 
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in the presence of quintessence. The boost factor is still of the order of ~ 10^ yfrj^ "t,ioo- In our fiducial illustration, 
we actually obtained an increase of the neutralino relic abundance by a factor of 3, 000 with toiqo = 2.5 and 77$ = 1 
in good agreement with the bench mark value which has been derived here. 

The increase of fi^ h? with 77$ has interesting consequences and brings up new perspectives as regards neutralino 
dark matter. To commence, the various avatars of the minimal or non-minimal supersymmetric extensions of the 
standard model start to be constrained - should R parity be conserved - by the accelerator data on the one hand 
side and by the requirement that the neutralino relic abundance should not overdose the universe or even exceed the 
observed value of — 0.135 ± 0.009 [15]. If a period of kination takes place in the pre-BBN period, the various 

SUSY configurations in the {^^h? ^m^) plane that arc so far allowed are shifted upwards with the consequence of 
becoming forbidden. Exploring in greater detail this question is a worthwhile project. 

We already anticipate that configurations with a very small relic density - for instance those for which poles dominate 
in the annihilation mechanism - would become cosmologically attractive if 77$ is large enough. The difference with 
the conventional cosmology lies in the significant enhancement of the annihilation cross section of neutralino dark 
matter candidates. At fixed Q,^ h? , notice that CTan^^ increases precisely by the same factor of ~ 10^ 'tiioo which 
we have derived above. This means a general enhancement of the various indirect signatures for supersymmetric 
dark matter. If neutralinos dominate the mass budget of the Milky Way halo, they should still annihilate today and 
produce gamma-rays, antiprotons and positrons which may be detected through the corresponding distortions in the 
various energy spectra. As a matter of fact, the recent HEAT experiment has confirmed [16] an excess around 8 GeV 
in the positron spectrum of cosmic rays. A large boost factor of ~ 10^ - 10^ in most of the supersymmetric parameter 
space is needed to explain that excess in terms of a homogeneous distribution of annihilating galactic neutralinos. 
A certain degree of dumpiness is actually expected in most of the numerical simulations but even in the extreme 
case of [17], it does not exceed a few hundreds. A period of kination in the early universe could provide an alternate 
explanation for that boost factor. 

Another potential consequence of quintessence is the rehabilitation of a fourth generation heavy neutrino in the 
realm of the dark matter candidates. In the conventional cosmology, a 100 GeV stable neutrino provides today a 
contribution of ^ lO""* to the closure density. Once again, kination at the time of decoupling would enhance that 
relic density and make that species cosmologically relevant. 

Notice finally that scenarios with extra-dimensions have the same effect as kination. The expansion rate is also 
increased and may even evolve as - to be compared to a behaviour in our case and to a dependence 
in the conventional radiation dominated universe. Implications of such scenarios on neutralino dark matter should 
be investigated. In the case of a low reheating temperature at the end of inflation, neutralinos are not thermally 
produced. Depending on the details of the scenario ~ decay of an inflaton field [18] or on the contrary decay of moduli 
fields [19] - the relic abundance is decreased or increased. 

A key ingredient of our study is the contribution 77$ of the quintessential scalar field to the overall energy density at 
the onset of BBN. A detailed analysis of the light element yields in the presence of kination [20] is mandatory at that 
stage in order to explore a promising scenario or to derive constraints on 77$ . The existence of dark energy opens up 
an exciting line of research and the study of its implications on the astronomical dark matter problem will certainly 
bring surprising results. 
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